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The energy spectrum of superfluid turbulence in the absence of the normal 
fluid is studied numerically. In order to discuss the statistical properties, we 
calculated the energy spectra of the 3D velocity field induced by dilute and 
dense vortex tangles respectively, whose dynamics is calculated by the Biot- 
Savart law. In the case of a dense tangle, the slope of the energy spectrum is 
changed at k = 2Tr/l, where I is the intervortex spacing. For k > 2tt/1, the 
energy spectrum has behavior in the same manner as the dilute vortex 
tangle, while otherwise the slope of the energy spectrum deviates from k~^ 
behavior. We compare the behavior for k < 2tt/1 with the Kolmogorov law. 

PACS numbers: 67.40.Db, 67.40. Vs. 

1. INTRODUCTION 

The particular attention is recently focused on the similarity between 
superfluid turbulence and conventional turbulencellili. In conventional tur- 
bulence, the energy spectrum E{k) has the Kolmogorov form k~^/^ for 
wavenumbers k in an inertial range. Such spectrum was observed indirectly 
by Stalp et al. in superfluid ^He above 1.4 in spite of the fact that the 
rotational flow in the superfluid component must take the form of discrete 
quantized vortex lines, and that there can be no conventional viscous dis- 
sipation in the superfluid component. This is understood by the idea that 
the superfluid and the normal fluid are likely to be counled together by the 
mutual friction and to behave like a conventional fluidEla. Since the density 
of the normal fluid is negligible at mK temperatures, an important question 
now arises: even free from the normal fluid, is the superfluid turbulence still 
similar to the conventional turbulence or not? 

Recently Nore et al. solved the 3D Gross-Pitaevskii equation and dis- 
cussed similarity to the Kolmogorov lawu. However this simulation includes 
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such complicated compressible effects as the radiation of sound from the 
vortex lines, the interaction between vortex lines and sound, etc. In order 
to consider pure superfluid turbulence in a simpler situation, we study the 
energy spectrum of the 3D velocity field induced by the vortex tangle under 
the vortex filament model which describes the dynamics of incompressible 
fluid. 



2. NUMERICAL CALCULATION 

For superfluid ^He, the vortex filament model is very useful, because 
the vortex core radius a ~ 10~^cm is microscopic and the circulation k = 
9.97 X lO^^cm^/sec is fixed by quantum constraint. Thus the dynamics of 
vortices is calculated numerically by the Biot-Savart law which is described in 
our previous paperi. In our calculation, a vortex filament is represented by a 
single string of points at a distance apart. When two vortices approach 
within A^, it is assumed that they are reconnected. The computational 
sample is taken to be a cube of size L =lcm. This calculation assumes the 
walls to be smooth and takes account of image vortices so that the boundary 
condition may be satisfied. 

We will introduce the energy spectrum of the velocity field induced by 
the vortex filament. Using Parseval's theorem / dk\v(k)\'^ = (27r)~^/ dr\v{r)\'^, 
the kinetic energy can be written as 

E = ^J dr\v^{r)\' = f (2vr)3/ dk\v{k)\\ (1) 

where ps is the superfluid density. By using the relation v{k) = ik x 
a>(fc)/|fcp, we obtain the kinetic energy: 

where 4j( A;) is the Fourier component of the vorticity a;(r) = Vx7;(r). In the 
vortex filament model, the vorticity can be defined as uj{r) = k J d^s'{$,)d{s{^) — 
r), so that a>(fc) is given by 

'^^^^ = j£rl dCe-''^^^"s'{0- (3) 

A vortex filament is represented by the parametric form s = where 
s refers to a point on the filament, the prime denotes differentiation with 
respect to the arc length ^ and the integration is taken along the filament. 
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The energy spectrum E{k) is defined as E = dkE{k). Finally we obtain 
the energy spectrum: 

E{k) = ^1^11 d^id^^AC,) . .'(6)e-^-(^(^^)-«^)), (4) 

where dQk denotes the volume element k'^ sinOkdOkdcjik in the spherical coor- 
dinate. The energy spectrum E{k) is calculated for the vortex configuration 
s(^) obtained by the simulation of the dynamics. 



3. ENERGY SPECTRUM OF VORTEX TANGLE 

In order to make clear the statistical properties of superfluid turbulence, 
we calculated the energy spectra of the 3D velocity field for dilute and dense 
vortex tangles respectively. 

First we discuss the energy spectrum for a dilute vortex tangle. This 
calculation of the dynamics is made by the space resolution = 4.58 x 10~^ 
cm and the time resolution At = 6.25 x 10~^ sec. Figure 1(a) shows the 
initial configuration of four vortices. Four rings move toward the center of 
the cube by their self-induced velocity to make the reconnection. After the 
reconnection the four rings move outside oppositely (Fig. 1(b), (c), (d)). 
The time evolution of the energy spectrum is shown in Fig. 2. The slope of 
the spectra is not strongly affected by the vortex configuration. The energy 
spectrum E{k) is prooortional to k^^, which comes from the velocity field 
of an isolated vortexEJa. 

Next the energy spectrum for a dense vortex tangle is discussed. This 
calculation of the dynamics is made by the resolution A,^ = 1.83 x 10^^ and 
At = 1.0 X 10^'^sec. Figure 3(a) shows the initial configuration of the vortex 
tangle where the direction of the circulation of upper vortices is opposite to 
that of lower vortices(Taylor-Green vortexi). The vortices become tangled 
through lots of reconnections(Fig. 3(b), (c) and (d)). The energy spectra 
are shown in Fig. 4. In Fig. 4(a) the spectrum may reflect an artifact of the 
initial configuration. However through the chaotic dynamics, these spectra 
reach the equilibration(Fig. 4(d)). Figure 5 shows that the time averaged 
spectrum of five conflgurations around t = 73.0 sec. The slope is changed 
at /c ~ 20 cm~^ and at k = 2tt/1, where / is the intervortex spacing. For 
k > Itt/I, the energy spectrum has k~^ behavior in the same manner as the 
dilute vortex tangle, by the contribution of the isolated vortex lines. For 
2'k/L < k < 20 cm~^, the spectrum is proportional to k~^^'^ by the effect 
of boundary. The spectrum for 20 cm^^ < k < 27r// is associated with 
velocity fleld at scales larger than the intervortex spacing, and this slope is 
very similar to the Kolmogorov law. 
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Fig. 2. Evolution of the energy spectrum of four vortices: t=0 sec(solid), 
t=6.25 sec (long-dashed), t=12.50 sec(dot-dashed) and i=18.75 sec(dashed). 
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(c) (d) 
Fig. 4. The energy spectra of the vortex tangle at t=0 sec(a), i=30.0 sec(b), 
i=50.0 sec(c) and i=73.0 sec(d). 
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Fig. 5. The time averaged spectrum of five configurations around t = 73.0 
sec. 



4. CONCLUSIONS 

This work studies numericahy the energy spectrum of the superfluid 
turbulence without the mutual friction. For k > 2Tr/l, the spectrum can be 
attributed to the contribution of the individual vortex lines. For 2tt/L < 
k < 20 cm~^, the spectrum depends on the boundary. In the intermediate 
range 20 cm~-^ < k < 27r/l, the slope of the spectrum is very similar to the 
Kolmogorov law. The calculation of the spectrum of a denser tangle is in 
progress so that the intermediate range may be extended. 
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